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Abstract. This paper investigates different Poisson structures that have been proposed 
to give a Hamiltonian formulation to evolution equations issued from fluid mechanics. Our 
aim is to explore the main brackets which have been proposed and to discuss the difficulties 
which arise when one tries to give a rigorous meaning to these brackets. Our main interest 
is in the deflnition of a valid and usable bracket to study rotational fluid flows with a free 
boundary. We discuss some results which have emerged in the literature to solve some of the 
difficulties that arise. It appears to the author that the main problems are still open. 



1. Introduction 

The aim of this paper is to present a discussion of numerous attempts to use the Hamil- 
tonian formalism of classical mechanics in hydromechanics and especially in the study of 
water waves. My motivation for this critical review came after a common work with David 
Sattinger [16] and some discussions with Adrian Constantin about some of his work on water 
waves with vorticity [10, 9]. 

The interest for this subject goes back to Zakharov [27] who showed that irrotational 
gravity waves could be given a Hamiltonian canonical structure. It was also influenced by 
the success of the Hamiltonian formulation for one dimensional evolution equations such as 
the Korteweg-de Vries equation, a theme which has been extremely intensive in the seventies. 

The Hamiltonian structure we refer to in this paper is that of general Poisson brackets 
which gives a more general framework in the sense that Hamiltonian systems can be deflned 
which are not necessarily canonical. If this structure is well understood on flnite dimensional 
manifolds, it is not the case for functional spaces. These structures have been deflned at a 
formal level, in the context of variational calculus [12, 13]. In inflnite dimension, the brackets 
are not deflned for all "smooth functionals" as it is the case for Poisson brackets on flnite 
dimensional manifolds, but only for a subclass of such functionals. This leads therefore to 
two natural questions: is the bracket closed for the class of functionals on which it is deflned 
and is the Jacobi identity^ satisfled by this bracket ? 

It appears that until recently, these questions have not been considered carefully and that 
Poisson brackets in functional spaces were deflned up to boundary terms as has been pointed 
out by Soloviev [23] for instance. 

It is however a fundamental question to check that a proposed bracket is a valid Hamil- 
tonian structure if one intends to go further than just a formal rewriting of the equations. 
To illustrate this fact I will just quote the work of Arnold [1] who was able to formulate a 
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^The Jacobi identity is the fundamental equation which must be fulfilled by a Poisson bracket: 
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stability theorem for plane flows using a method now known as the Energy- Casimir method. 
This work relies on the existence of Casimir functions^ for the underlying structure. 

Besides, one could suggest that the ultimate goal of introducing Hamiltonian formalism 
in hydrodynamical problems and especially in the study of water waves would be to derive 
from it new results (like e.g. Arnold's stability theorem [1] or some recent results on particle 
trajectories [5, 6]). To achieve this, one cannot however avoid the difficult question of defining 
a valid Hamiltonian structure. 

This paper proposes to discuss this question with a critical review of the main Hamiltonian 
structures which have been proposed in the literature (up to the author's knowledge). 

The content of the paper is as follows. In Section 2, we review the basic material on 
Poisson structures for finite dimensional manifolds. In Section 3, we extend these definitions 
to functional spaces and raise the main difficulties which appear when one tries to define 
valid brackets in this more general context. In Section 4, we discuss Hamiltonian structures 
on the space of smooth functions on the circle, where things work well. Section 5 is devoted 
to Arnold's bracket, a formulation of the Lie-Poisson bracket for the Lie algebra of divergence 
free vector fields on a compact domain and which is the background structure for the motions 
of an ideal fiuid with a fixed boundary. Several versions of this bracket are proposed and 
discussed. In the final section. Section 6, we discuss some brackets which were introduced in 
[17] to study the difficult problem of fiuids with vorticity and free boundary. It appears that 
the proposed bracket is not closed. 



2. Poisson brackets in finite dimension 

2.1. Symplectic and Poisson manifolds. A symplectic manifold is a pair (M, tu), where 
M is a smooth manifold and u; is a closed and nondegenerate 2-form on M. Such structures 
appear naturally in mechanics (see [14]). If is the configuration manifold of a mechanical 
system, its phase space is the cotangent bundle T*N and is equipped with the canonical 
2-form given by: 

"^dpi A dq\ 

i 

Since a symplectic form u is nondegenerate, it induces an isomorphism TM T*M. The 
inverse of this isomorphism defines a skew-symmetric bilinear form P on the cotangent space 
T*M and a skew-symmetric bilinear mapping on C°°(M), the space of smooth functions 
/ : M ^ M, given by 

(1) {f,g} = P{df,dg), f,geC°^{M), 

called the Poisson bracket of the functions / and g. For example, when M = T*N is a 
cotangent bundle, the corresponding bracket, known as the canonical bracket is given by: 

If ^ 9f dg dg df 

Z-j Qp. Qqi Qp. Qqi 

The observation that a bracket like (1) could be introduced on C°°(M) for a smooth 
manifold M, without the use of a symplectic form, leads to the general notion of a Poisson 
structure. 



Casimir function is a smooth function whose bracket with every over smooth function vanishes. 
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Definition 2.1. A Poisson structure on a smooth manifold M is a skew-symmetric bilinear 
mapping {f,g) ^ {f,9} on the space C°°(M), which satisfies the Jacobi identity 



Each Poisson bracket {, } corresponds to a smooth field P of bivectors, called the Poisson 
bivector of (M, {, }) and such that 



for all f,g E C°°{M). The Jacobi identity implies that the bivector field P must satisfy a 
certain condition, namely that [P, P] = 0, where [ , ] is the Schouten-Nijenhuis bracket"^. 
The Hamiltonian vector field of a smooth function / on M is defined by 



These functions play an important role in the study of the stability of equilibrium of Hamil- 
tonian vector fields. Notice that in the symplectic case, the only Casimir functions are the 
constants. 

2.2. Poisson reduction. Let us now explain how these Poisson structures appear naturally 
in mechanics. Let the configuration manifold of a mechanical system and M = T*N its 
corresponding phase space. It often happens that the system has some symmetries repre- 
sented by the (left) action of a Lie group G on A^. This action lifts to a symplectic action 
of G on M = T*N, that is each diffeomorphism induced by an element g G G is a canonical 
transformation of M = T*N. If the group G acts freely and properly on M, the reduced 
phase space M/ G is a manifold and we may ask which structure from M is inherited by the 
quotient space M/G. 

For that purpose, let tt : M M/G be the canonical projection. Notice that ker7r'(a:) is 
the tangent space to the G-orbit through x. Let u; be a 2-form on M, P a bivector field on 
M and recall the following criterions 

(1) There exists a bivector field P on M/G such that n' o P = P o n if and only if 




{{f,g},h} + {{g,h},f} + {{h,f},g} = 0, 



{f,gh} = {f,g}h + g{f,h}. 



{f,g} = P{df,dg) 



Xf = Pdf 

so that {f,h} = Xh ■ f. The Jacobi condition on P insures that 

as in the symplectic case. 

A Casimir function is a smooth function G on M such that 

{G,/} = 0, V/gG~(M). 



{g*P){x) - P{x) G ker7r'(a;) 

for each point x G M. 
(2) There exists a 2-form u on M/G such that h^lj = u if and only if 

g*LJ = u and ix^jJ = 

for each vector X G keiTr'. 



The Schouten-Nijenhuis bracket is an extension of the Lie bracket of vector fields to skew-symmetric 
multivector fields, see [26]. 
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Notice that, unless G is a discrete group, the second condition on uj is never satisfied and 
hence the symplectic structure on M = T*N cannot get down to M/G. However, condition 
(1) is fulfilled by the Poisson bivector P of any Poisson structure on M invariant under 
G, and leads naturally to the existence of a reduced Poisson structure on M/G such that 
vr : M — i> M/G is a Poisson map, i.e. such that 

{/o7r,5(0 7r} = {/,5-}o7r 

for all f,g & G°°{M/G). This process is known as the Poisson reduction [18]. 

2.2.1. Lie-Poisson structure. The main illustration of this reduction process leads to the Lie- 
Poisson bracket. Let G be a Lie group and Q its Lie algebra. The left action on G lift to a 
symplectic action on T*G ~ G x g* (equipped with the canonical symplectic structure) and 
induces a Poisson structure on T*G/G ~ g* given by 

(4) {/, g}im) = -m{[dmf, dmg]) 

for m G g* and f,g E G°°(g*)''. The corresponding Poisson bivector P is given by 

Pra{df,dg) = ad*^fm{dg) 

where ad* is the coadjoint action of g on g*. 

2.2.2. Euler equation. The Lie-Poisson structure is the framework for the evolution equa- 
tion known as the Euler equation on a Lie group G. Consider a one-sided (left or right) 
invariant Riemannian metric < ■, ■ > on G. The geodesic fiow corresponds to the flow of the 
Hamiltonian vector field on T*G equipped with the canonical structure and Hamiltonian 

The reduced Hamiltonian function Ha and the reduced Hamiltonian vector field Xa on g* 
are 

HA{m) = - (m, A~^m), XA{m) = ad\-i^m, meg* 
where A : u h-^< u, ■ >e is called the inertia operator. 

Example 2.2 (The rigid body). Euler equations of motion of a rigid body: 
h -h h- 1 1 h - h 

iOl = UJ2UJ3, U2 = tUl^S, ^^3 = } t^lt^2 

^2 ^3 

are the basic example of Euler equations. In that case, the group G is the rotation group 
SO (3). The Lie-Poisson bracket on so (3)* ~ is given by 

{/,^}M = m- (grad/(m) Agrad^(m)), f,ge G'^{R^), 

and the Hamiltonian is 

where h^h, I3 are the principal moments of inertia of the rigid body and rrik = Ik^k- 

^Here, c?m/, the differential of a function / G C°°(g*) at to G g* is to be understood as an element of the 
Lie algebra g. 
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3. POISSON BRACKETS IN FUNCTIONAL SPACES 

Several authors have tried to extend the notion of Poisson brackets to functional spaces in 
order to study evolution equations, see [20] for an excellent overview of the subject. There are 
however serious difficulties to handle when one enters into the details of these constructions 
as was pointed out in [21, 22, 23, 24, 25]. 

In this section, we will review some of these difficulties. We consider Poisson brackets 
for smooth functionals defined on the the space C°°{M) of smooth functions on a manifold 
M or more generally on the space of smooth sections T{E) of a vector bundle over M (for 
simplicity, we will suppose that M is the closure of an open subset of the Euclidean space 
M"' with smooth boundary). 

3.1. Directional derivative versus variational derivative. Let F be a smooth real func- 
tion on some Frechet vector space C'^{M,E) where -E is a finite dimensional vector space. 
The directional derivative or Frechet derivative of F at u in the direction X G C°°{M^E) is 
defined as 

d 



DxF{u) = - 
as 



F{u + eX). 

£ = 

In general, the directional derivative X ^ DxF{u) of a smooth functional F is nothing more 
than a continuous linear functional on C°°(M, F). Sometimes, this linear functional can be 
represented as 



f 6F 

DxF{u)= / —{u)-XdV, VXeC°°(M,F) 
Jm 



where 



5F 
ou 



is a smooth map (vector field) from C°°(M, F) to C°°(M, F). The vector field 6F/6u is 
unique and we call it the gradient of F. 

There is another notion of derivative, whose origin comes from variational calculus 

d 



DF(u).Su 

de 



F{u + e5u) 



=0 



where the variation 6u has compact support and is subject to various boundary conditions. 
We call it the variational derivative of F. At first, it seems that the two definitions are the 
same. Of course, this is the case if M is a compact manifold without boundary, but in general 
it is not. 

A function F on C°°{M,E) is called a local functional if 



F{u) = / f{x,u^''^)dV 
Jm 

depends of u through a smooth function / (the Lagrangian density of F) which depends only 
on X and the r-jet of u up to a certain order r. In that case, the Frechet derivative of F is 

F(u + e6u)= [ y"^(x,u'''-^)6u^if\x)dV 



DsuFiu) = ^ 
de 



^=0 J,fc 

where ui, . . . ,Up are the components of u and 
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Using the Leibnitz rule repeatedly [20], we can show that 




where is the Euler operator defined by 

J ^ 
P is a (functional) vector field 

P(x,mW) = (Pi(x,mW),...,P„(x,m(^))) 

and the divergence of P is defined by 

div P = DiPi + D2P2 + ■■■ + DnPn, 

where Di = d/dx'' is the total derivative with respect to x\ 

Therefore, the variational derivative of a local functional F can always be put in a gradient 
form 

DF{u).Su= 1 SF-SudV 

J M 

where 

5P=(Ei(/),...,E,(/)). 

However, when the manifold M has non-empty boundary, the variational derivative and 
the Frechet derivative may differ by a boundary term. A local functional does not have 
necessarily a gradient relatively to its Frechet derivative. 

Example 3.1. This may happen for instance for a local functional given by 

F{u)= / div PdV= / P-ndS 

J M JdM 

The variational derivative of F is identically zero but the Frechet derivative of F has no 
reason to vanish and cannot be put into gradient form. This problem arises because in 
the definition of the Frechet derivative we allow all smooth variations whereas in the definition 
of the variational derivative we allow only variations subject to boundary conditions. 

A Poisson bracket {F, G } is first of all a bilinear map depending on the first derivative of 
F and G. Contrary to the finite dimensional case, it seems extremely difficult to define a 
tractable Poisson bracket on the set of all functionals. The reasonable thing is to restrict the 
definition of the bracket to a subclass of functionals. For instance, in the formal variational 
calculus [20], a Poisson bracket is defined on the subclass of local functionals through a 
bilinear map on their variational derivatives but this bracket appears to be defined up to 
divergence terms. When the manifold is compact without boundary this may lead to a 
coherent Poisson bracket but when the manifold has non-empty boundary some difficulties 
arise. 

Example 3.2 (The Gardner bracket). It was discovered by Gardner, [11], that the Korteweg-de 
Vries equation 

Uj^ Uxxx ~t~ ^^x 
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can be written as a Hamiltonian equation using the bracket 



{F,G}iu)= [ ^J^nJ-^dx. 

Jgi ou ou 



and the Hamiltonian 



H{u) = [ f-^w' + ^w' J dx. 



3.2. Closure of the Poisson bracket and Jacobi identity. As we have just seen, there is 
no well-defined Poisson bracket on the space of all smooth functionals. The known brackets 
are defined on a subclass A of functionals, called admissible functionals. 

When the manifold M is compact without boundary, it is possible to choose for A the 
whole space of local functionals. We may then define a Poisson bracket {F, G } on ^ using 
an expression like 



where P is a linear differential operator (witch may depend of the r-jet of u), as in the 
Gardner bracket. This gives us a well-defined bilinear map 

Ax A-^ A 

since the expression we have for {F, G } is itself a local functional. 

When M has non-empty boundary this is not sufficient and some other boundary conditions 
have to be introduced (see Section 5). Now this leads to an immediate other question: If F 
and G satisfy this boundary conditions, is this true for {F, G}? In other words is the class 
A of admissible functionals (verifying the boundary conditions) closed under the bracket ? 
As we shall see this is not at all obvious. 

Finally and last but not least, if all these required conditions are satisfied, we have to check 
that the bracket verifies the Jacobi identity 



This last verification can be very tedious but the real difficulty remains however the closure 
of the bracket. 

3.3. Hamiltonian structures. All these considerations lead us to introduce the following 
scheme to define a Poisson bracket on a functional space. First define a subspace A of 
smooth functionals (local functionals for instance, if dM = 0, or local functionals with some 
boundary conditions otherwise). Then we introduce the following definition of a Hamiltonian 
structure^ on A. 

Definition 3.3. A Hamiltonian structure on ^ is a bilinear operation {-, ■ } on ^ such that 
for any F,G, H E A we have: 



Remark 3.4. Notice that the Leibnitz rule has been eliminated from the definition of a Hamil- 
tonian structure. In fact, there is no well-defined commutative product on local functionals. 

^The terminology Hamiltonian structure is commonly used instead of Poisson structure for functional 
spaces. 




{{F,G},H} + {{G,H},F} + {{H,F},G} = 0. 



(1) {F,G}e A 

(2) {G,F} = -{F,G}, 

(3) {{F,G},H} + {{G,H},F} + {{H,F},G} = 0. 
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In the following sections, we review some well-known brackets that have been proposed in 
the literature (see also [15, 8]). 



4. The Lie-Poisson bracket on Vect*(S'^) 

In this section we will consider the Lie-Poisson bracket on the "dual" of the Lie algebra of 
smooth vector fields on the circle Vect(S'^) ^ C°°(S'^). Recall that the canonical Lie-Poisson 
structure on the dual q* of a Lie algebra g is given by 

{F,G}{m) = -m{[dmF,dmG]). 

To give a sense to this expression, we have first to define an injection from q to q*. 



4.1. The regular duaL Since the topological dual of the Frechet space Vect(S'^) is too big 
and not tractable for our purpose, being isomorphic to the space of distributions on the circle, 
we restrict our attention in the following to the regular dual g*, the subspace of Vect*(S'^) 
defined by linear functionals of the form 

Ml—* / mu dx, 

for some function m G C°°(5'^). The regular dual q* is therefore isomorphic to C°°{S^) by 
means of the inner product^ 

< u,v >= / uv dx. 



4.2. Local functionals. A local functional F on Vect*(S'^) ~ C°°{S^) is given by 

F{m) = I f{x,m,mx,...,m^^^)dx. 

Since there are no boundary terms, its functional derivative DF{m) is equal to its variational 
derivative 

r 

DFim).5m = / —6mdx, meC^(S^). 
J si dm 

where 

- = Ti-D.y^. 

j=0 

The map m SF/Sm can be considered as a vector field on C°°{S^), called the gradient of 
F for the L^-metric. In other words, a local functional on C°°(5'^) has a smooth gradient. 



"In the sequel, we use the notation u,v, . . . for elements of g and m,n, . . . for elements of g* to distinguish 
them, although they all belong to C°°{S^). 
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4.3. Hamiltonian structures on the regular dual. To define a Poisson bracket on the 
space of local functions on Vect*(5'^), we consider a one-parameter family of linear operators 
Pm ijn G C°°(5'^)) whose coefficients are smooth function of x, m and a finite number of its 
derivatives and set 

(5) {F,G}(m)= f 5FP^5Gdx. 

where SF and SG stand here for the variational derivatives 6F/Sm and SG/Sm. The operators 
Pm must satisfy certain conditions in order for (5) to be a valid Hamiltonian structure on the 
set A of local functionals on the regular dual Vect*(S'^). First it must be a skew-symmetric 
operator (relatively to the inner product). 

/ SFPm5Gdx = - [ 5GPmSFdx, \/F,GeA. 

Since the expression for {F, G } is a local functional, the class of local functional is closed 
under this bilinear operation. Therefore we need only a criteria on P to ensure that Jacobi 
identity is satisfied, in order to obtain a Hamiltonian structure. 

Lemma 4.1. The Jacobi identity for (5) is equivalent to the condition 

(6) 0/ SF {DpsHP)SGdx = 

for all F^G^ H G A where O indicates the sum over cyclic permutations of F, G, H and DsmP 
is the Frechet derivative of P in the direction 6m. 

Remark 4.2. Notice first that since P is a linear differential operator whose coefficients are 
smooth functions of a;, m, m^,, . . . , the Frechet derivative of P in the direction 6m is just 
the linear differential operator obtained by replacing the coefficients of P by their Frechet 
derivatives in the direction X. Since P is assumed to be skew-symmetric, so is DsmP- 

Proof. We already know that {F, G } is a local functional and hence its variational derivative 
6{F, G } is an gradient for {F, G }, that is 

Dsm{F,G}= [ 6{F,G}6mdx. 

By definition of the bracket, we have 

{{F,G},H}= [ 6{F,G}P6Hdx = Dpsh{F,G}. 

Using the fact that the second Frechet derivative is a symmetric operator and the fact that 
P is skew-symmetric, we get 

Dsm{F, G}= [ [{Dsm6F) P6G - {Dpsf6G) 6m + 6F (DsmP) 6G] dx 

and hence 

{{F,G},H}= [ [{Dpsh6F) P6G - {Dpsf6G) P6H + 6F (DpshP) 6G] dx. 

Taking the sum over cyclic permutations of F,G,H, the two first terms of the right hand 
side of the last equation cancel and we obtain the equivalence of Jacobi identity with (6), 
which ends the proof. □ 
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To check (6) is still tedious in practice. Following Olver [20], it is preferable to use the tech- 
nique of functional bivectors, which generalizes the notion of Poisson bivectors and Schouten- 
Nijenhuis brackets. First, given a functional density /(x, m^'')), define 

where stands for the total derivative relative to x. Extending the action of the differential 
operator P on 6' in an obvious way, we can write 

{F,G} = ]-f {9{5F){P9){5G)-9{5G){P9){5F)} dx=l [ {9 A P9){5F,5G) dx 

so that 

Q = - [ {9AP9}dx 
2 Jsi 

appears as the analogue of the Poisson bivector for finite dimensional Poisson brackets. 
Example 4.3. For the Gardner bracket we have 

Q = l [ {9 A 9,,} dx. 

Proposition 4.4 (Olver [20]). A skew- symmetric linear differential operators P (with coef- 
ficients depending on x, m, m^, ...) defines a Hamiltonian structure on the space A of local 
functionals on Vect*(S'^) if and only if it satisfies 



[ {9 A{DpeP A9)}dx = 0. 



Remark 4.5. Notice that the preceding expression is an alternatinng trilinear expression on 
functional densities. Note also that the two wedges have different meanings. The first one 
corresponds to wedging the ordinary multiplication of two functional densities whereas the 
second one is the wedging relative to the non-abelian bilinear operation {f,g) ^ {DpfP)g. 

Proof. Let F, G, H be local functionals and 6F, 6G and 6H their variational derivatives. Then 

- / A {DpeP A 9)} {5F, 5G, 5H)dx = [ SF (DpshP) 5G dx. 
2 Js^ Js^ 

Hence the proposition is just a corollary of lemma 4.1. □ 

Example 4.6. The Gardner bracket or more generally the bracket obtained from a skew- 
symmetric differential operator P with constant coefficients satisfies the Jacobi identity since 
the Frechet derivative of such operators in any direction is zero and hence DpgP = 0. 

Example 4.7. The canonical Lie-Poisson structure on Vect*(5'^) is given by 

(7) {F, G}{m)= I m [5F, 6G] = - [ SF {mD + Dm) SG dx 

It is represented by the skew-symmetric operator 

P = - (mD + Dm) = - {2mD + mj) 

where D = d/dx. We get 

DpeP = {Am9^ + 2m^9) D + {2m9^^ + 2m^9^ + m^^9) L 

hence 

DpgP A9 = 2m J A 9^ + 2m9^^ A9 + 3m A 9 
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and 



9 A {DpgP A9) = 0. 
Example 4.8 (Burgers equation). The inviscid Burgers equation 

Ut = —3uUx 



can be written as an Euler equation on Vect*(5'^) with the Lie-Poisson bracket (7). It 
corresponds to the inertia operator m = Au = u and Hamiltonian 



Example 4.9 (Camassa-Holm equation). The Camassa-Holm equation [4] 

Ut — Utxx + 3nn^ — 2uxUxx — uUxxx = 

can be written as an Euler equation on Vect*(5'^) with the Lie-Poisson bracket (7). It 
corresponds to the inertia operator m = Au = u — Uxx and Hamiltonian 



cf. [19] - see also the discussion in [7]. 

Notice however that H is not a local functional of m since it depends on m by the inter- 
mediary of the non local operator A~^. To overcome this difficulty, one may try to extend 
the Hamiltonian structure (7) for functionals which are local expressions x,u,Ux,m,mx, . . . 
where u = A~^m rather than x, m, nix, ■ ■ ■ ■ But this space of functionals is not closed under 
the preceding bracket and the space of functionals generated by successive brackets of such 
functionals seems tedious to describe. 

In that case however, it is possible to overcome these difficulties by extending the Hamil- 
tonian structure to the whole space of smooth functionals which have a smooth gradient, 
that is 



where m i— > 5F{m) is a smooth smooth map from C°°{S^) to C°°{S^). The first two axioms 
which define a Hamiltonian structure are verified. Indeed, the bracket of two such functionals 
has itself a smooth gradient, namely 



Finally, Jacobi identity is also verified for this extension. In fact. Lemma 4.1 and Proposi- 
tion 4.4 are still valid for those more general functionals. 



Let n be a relatively compact domain in or with a smooth boundary. We let 
SDiff(fi) be the group of volume-preserving smooth difi^eomorphisms of ^2 and SVect(fi) the 
Lie algebra of divergence-free vector fields on f2, tangent to the boundary, which can be 
interpreted as the Lie algebra of SDiff(f2). In a famous article [2], Arnold showed that the 
Euler equations of perfect incompressible fiuid fiows in a fixed domain 

du 

(8) — — h VuM = — gradp, divw = 0, m ■ n = on dVL 






6{F, G } = DpsF SG - DpsG SF + SG DxSF - 6F DxSG. 
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could be interpreted as a the Euler equation of the right-invariant (weak) Riemannian metric 

(9) <u^,v^>= / u^-v^dV 



where u^,v^ are vector fields over ip G SDiff(r2) {Lagrangian velocities). 

The regular dual of SVect(r2), noted SVect*(fi), consists of linear functionals on SVect(r2) 
with smooth density a G r2^(fi) 

a(u) dV. 



Since exact one-forms are L^-orthogonal to divergence-free vector fields, a one-form a G f2^(f2) 
represents an element of SVect*(f2) only up to total differential. Each element of SVect*(r2) 
is therefore represented by a class [a] in fl^ (fl) / dfl^ (fl) . If moreover, is a simply connected 
domain, (Q) / dfl^ (Q) is isomorphic to dfl^{fl) via the exterior derivative d and the class 
[a] is completely represented by the two-form u = da, called the vorticity. 

5.1. Arnold bracket. The bracket, now known as Arnold bracket, is defined for smooth 
functionals F on SVect*(r2) whose Frechet derivative can be written as 

(10) D[Sa]F{io) = 5a (^^^ dV where ^ G SVect(l]), 

in other words, for those functionals which have a gradient in the Lie algebra SVect(n). 
It is given by the formula^ 

^6£ 6G' 

Suj ' 6uj 



(11) {F,G}ico) = - / a 



M 



M 



\5uj 5uj ) 



Proposition 5.1. The bracket defined by equation (11) is a valid Hamiltonian structure on 
the space of smooth functionals on SVect*(f2) which have a smooth gradient in SVect(f2). 

Proof. We have to check the three properties of Definition 3.3. Expression (11) is clearly 
skew-symmetric. To show that the bracket is closed, we recall first that the symmetry of the 
second Frechet derivative leads to 



5a{Dysp]^F) dV= I 6f3{D[Sa]SF) dV 
Jn 



for every admissible functional F. From this property, we deduce that for every admissible 
function F and G, we have 

D[Sa]{F,G}{uj)= [ 6a{D[,,^^]6G-D[,,^^]6F-[SF,6G]) dV, 
Jn 

that is {F, G } is also admissible with gradient 

6{F, G } = D[i,^^]6G - D[,,^^]6F - [6F, 6G] . 
It remains to check Jacobi identity. We can write 

{F,G}{uj)= [ 5F{uj)- P^5G{u)dV 



^The equality of the two formulations results from 

da{u, v) = u ■ grada(u) — v ■ grada(w) — a{[u, v]). 
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where 

: SVect(fi) ^ SVect(n), u ^ -Proj (i„cj) 

and Proj is the projection from onto SVect(fi), which associates to a one- form a the 

unique divergence free vector field tangent to the boundary, such that 

I a{w) dV= [ v-wdV, Vw e SVect(fi). 
Jn Jn 

With these notations, we get as in Lemma 4.1 

O {{F,G},H} = - O [ [SF, 6G] ■ P5HdV 

= ! u{6H,[6F,6G])dV 



O / a{[6H,[6F,6G]])dV = 
Jn 



where u = da. □ 

5.1.1. Euler-Helmholtz equation. Arnold's bracket (11) can be used to interpret Euler's equa- 
tions of perfect incompressible fluid flows (8) in their Helmholtz or vorticity representation 

(12) dtuj = curl(M X u), u = curlw 

as the Euler equation of the metric (9). 

Recall that the curl of a vector held u is deflned as the unique vector u such that 

vol = dit 

where m'' is the covariant representation of u. Therefore, SVect*(f2), the space of exact two- 
forms can be identified with the space of curls and the inertia operator of the L? metric (9) 
can be described as 

A : SVect(fi) ^ SVect*(fi), u ^ curlw. 

This operator is invertible. Let u G SVect*(fi) be a curl. Then u = A"^uj is the unique 
solution of the problem 

curlw = uj, divu = 0, m ■ n = on dQ. 

The Hamiltonian is given by 

H{iu) = ^ ! WW dV. u = A-^u. 

We have 



D[Sa]H{u!) = u ■ 6udV, 6a = {6u)^ 
Jn 

and hence H is an admissible functional with gradient 

6H{uj) = u = A~^uj. 
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Euler equation F = {F,H}, for all admissible functional^ F gives 



dtu-5FdV= / iu- {6F xu)dV= / 5F-{uxu) dV, 
n Jq Jn 

that is 

dtu = u X u, modulo a gradient 

and taking the curl, we get 

dfUJ = curl(M X u). 

Remark 5.2. We could have restricted the definition of Arnold's Poisson bracket for local 
functionals which have a gradient. In fact this space is closed under the bracket. But this 
would not have permitted us to treat the hydrodynamic problem since the Hamiltonian is 
not a local functional (see Example 4.9 for a similar situation in dimension 1). 

Remark 5.3. In several papers, the Arnold bracket is written as 

{F,G}{u) = [ curln- x ^] dV. 



Jn \6u 5u J 

for smooth functionals with smooth gradient on the Lie algebra SVect(f2) rather than 
SVect*(f2). This is just the "pullback" of (11) by the inertia operator A. The fact that this 
bracket preserves the space of functionals which have a gradient is less obvious to see in 
this expression because of the term curlw which leads to an integration by parts, but in fact 
it works. The advantage of using variables u instead of u is that the Hamiltonian becomes a 
local functional in these variables. In that case, the Hamiltonian equation, F = {F, H } for 
all admissible F, leads directly to equations (8). 

Remark 5.4. A third interpretation of Arnold bracket was given in [20]. It is defined, in the 
context of formal variational calculus (where boundary terms are ignored) for local function- 
als on SVect*(fi). The gradient of a functional F is defined here as 

f 6F 

DF{uo).5uo= / —-dujdV 
Jn did 

where the gradient, SF/Suj is a divergence free vector field and the variation Su is assumed to 
vanish on the boundary. Notice that the definition of the gradient given here is quiet different 
from the previous definition (10). Indeed the two definitions differ through a boundary term 

5u; dV =[ Sa ( cuil^] dV + [ (^x5u\-ndS 



SuJ Jn \ J Jqq \ 6uj 

where 6uj = curl 6u and 6a = {6u)K Therefore we cannot conclude that both Poisson 
structure are rigorously equivalent. 

For two-dimensional fiows, the Hamiltonian operator P is represented as 

P = LOxDy — LUyDx 



^Each vector field u G SVect(f2) can be realized as the gradient of an admissible functional, namely of the 
linear functional 



F{uj) = / a{u) dV, 
Jq 



da 
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and the gradient of the Hamiltonian 



H{uj) = 11 \\u'^\\ dV, u = A^^uj. 

5H/ 5uj is the stream function ip of the velocity u. It was shown in [20] that in this context, the 
Jacobi identity was satisfied and that the Hamiltonian equation F = {F, H } was equivalent 
to Euler-Helmholtz equation (12). We insist on the fact that the computations which leads 
to these results relies on the vanishing of the variations on the boundary. 

5.1.2. Enstrophy. This Poisson bracket (11) has been rejected by the authors in [17] because 
for two-dimensional flows, the generalized enstrophy functional 

C{lj) = / (t){Lj) dx A dy 
Jn 

which is known to be invariant under the coadjoint action of SDiff(f2) on SVect*(f2) is not a 
Casimir function for this bracket. Indeed 

Disa]C{^) = / Sa[cml{(j)'{uj)k)) dx A dy + (p (f)'{uj)6a 
Jn Jan 

has some boundary terms and is therefore not an admissible functional for (11). 

5.2. Second LMMR bracket. Since Casimir functions play a fundamental role in the 
study of stability of two-dimensional flows as it has been shown in [1], the authors in [17] 
have proposed to improve the definition of Arnold's bracket by taking into account boundary 
terms so that the enstrophy becomes a Casimir function. 

They have derived this bracket using the same reduction process which has been used for 
Arnold's bracket. The difference lies in a different choice of admissible functionals. 

The starting point is the Lagrangian description of the problem. For an incompressible 
fluid moving in a flxed domain f2, the configuration space is the group of volume-preserving 
diffeomorphisms SDiff(r2). The phase space, T*SDiff(f2) has to be understood as the set of 
pairs (v?, /i) where ip G SDiff(r2) is the "base point" and /i is a one-form over ip (i.e. for each 

^, Ma^)er;(,)fi). 

The class of admissible functionals F, previously limited to smooth functionals which have 
smooth gradient, is now extended to ones whose "gradients" can be written as 

6F SF 

1^ — -J 1- — , -T— — -J h ^ao-j — , 

oif dip dip oji oji o/i 

where 5q^ is the Dirac measure on Vt concentrated on dVL ^. A Poisson bracket can be deflned 
for those functionals using the formal canonical bracket on T*SDiff(f2) 

{FG}= [ f^^-^^^ dV 

provided that the boundary condition 

(13) 6^6^ _ 6^6^ _^ 

6ip 6fi Sip Sfx 
is satisfled to avoid squares of delta functions. 



Notice however that this decomposition is not unique. 
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The Lie-Poisson reduction of the phase space T*SDiff(n) by the gauge group SDiff(J7) 
(corresponding to relabeUng fluid particles) leads to the Second LMMR bracket [17] deflned 
for functionals on Vect*(f2) whose Frechet derivative can be written as 

DF{u).5u = / -—■5udV+ / -—-SudS 
Jn JgQ du 

The expression for the resulting bracket is quite complicated and will not be given here. It 
must be stated, however that this bracket is well-deflned for a pair of admissible functionals 
[F, G) only if condition (13) is satisfled. This will be the case if one if one of the functionals F 
or G satisfy F/5u = 0. It was shown in [17] that for two-dimensional flows, the generalized 
enstrophy was a Casimir function for this bracket in the sense that {C, F} = for all 
functions admissible function F such that 5^ F/5u = 0. 

We will not try to check that this second LMMR bracket is a valid Hamiltonian structure. 
The definition of this bracket {F, G } requires the condition (13) on the pair of functionals 
(F, G) to be satisfied. But this latest condition concerns the pair (F, G) and not each of 
the functionals F, G alone. Therefore, it is not clear on which subclass of functionals is this 
bracket defined. 

5.3. Soloviev bracket. In a series of papers, [21, 24, 25], Soloviev tried to define a Poisson 
bracket for local functionals which avoids this tedious boundary condition (13). The idea 
introduced in [21] is to define a bracket involving not only the "first gradient" (the factor 
of Su) but the complete set of "higher order gradients" (the factor of (Su)^-^^) of a local 
functional. 

Using the Leibnitz rule but making no integration by parts, we can write the Frechet 
derivative of a local functional F as 



DF{u).6u= f V^(x,«M) 5ui'\x)dV= f TDj{Ei{f)6uk) dV 
J^ j,k ^""k Jn 

where the higher Eulerian operators^° Ej^ are defined by 



KDJ ^ ^ 

the binomial coefficients for multi-indices are 

K \ I ki \ ( k,. 



^ J ) \ Ji J \ J 
and 

(-d)^ = (-i)|^IDk. 

The following formula was derived by Soloviev in [21] to define a Poisson bracket on the class 
of all local functionals 

{F,G} = Y,Y. [ {EiU)Ip,E^{9)) dV 

J,K p,q 

where the operator Ipq are subject to certain conditions to satisfy Jacobi identity. 



^''Notice that all the sums are finite since only a finite number of derivatives appear in all these formula. 
The zero order higher Eulerian operator is just the ordinary Euler operator Ek. 
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Example 5.5. In [24], this method was illustrated for the formulation of Arnold's bracket 
presented in Remark 5.4 for 2 dimensional flows. The antisymmetric operator / was given 
in this case by 

1 

2' 

where the derivative of the characteristic function 9 = 9^ has to be understood in the sense of 
distributions using certain rules [24]. It was shown that, up to these rules, we obtain a valid 
Poisson structure. There is however one objection on this example: up to my understanding, 
Soloviev's formalism was developed for local functionals but the Hamiltonian giving rise to 
the Euler equations in this case is 



H{lj) = 2 ^ ''^"^ ^ ~ curlu, 



2 

||n|| 



which is not a local functional of the variable u. 

6. Poisson brackets for ideal fluids with a free boundary 

In 1968, Zakharov [27] showed that Euler's equations for irrotational gravity waves could 
be written as a canonical Hamiltonian system. The Hamiltonian is 

H = l JJJ{gTad^ydV + h JJe{x,y,t)dS. 

D M2 

The Poisson brackets implicit in Zakharov's observation are the canonical brackets 

the Hamiltonian flow is then the canonical flow 

The Hamiltonian H is regarded as a functional of () where ( = ({x,y,t) is the height 
of the free surface, and (p = (flan is the trace of the harmonic function (p on the free surface, 
with = on the bottom. The evolution takes place in the space of harmonic functions 
on D. Zakharov's result is verifled by calculating the gradients of H with respect to ( and 

In [17], a generalization of this Hamiltonian structure was proposed for incompressible fluid 
flows with possible vorticity. It is however no longer a canonical structure. The approach 
used in [17] to derive a Hamiltonian structure is essentially the same as the one used to 
derive Arnold's bracket (fixed domain): using a Poisson reduction process of the canonical 
symplectic structure on the phase space by the gauge group (relabelling of particles). The 
main difference is that in the free boundary case, the gauge group no longer acts transitively 
on the configuration space (the space of embeddings of a reference domain in M"). 

6.1. First LMMR bracket. This structure, known as the first LMMR bracket is defined 
on the space of pairs {v, S), where S is the free surface and v is the spatial velocity field, 
a divergence free vector field defined on Ds, the region bounded by S. The surface S is 
assumed to be compact and diffeomorphic to the boundary of a reference region D. 

The class A of functionals F : A/" — > M on which this bracket is defined is formed by 
functionals with the following properties: 
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(1) A variation 5v is just a divergence free vector field on Dj^ and we assume that there 
exists a divergence free vector field 5F/ 6v defined on D-^ such that for all variations 
Sv: 

r 6F 

D^F(v,J:) ■ 6v = / —-SydV 

where DyF is the derivative of F holding S fixed. 

(2) A variation which is a function on S, has to be understand as an infinitesimal 
variation of S in its normal direction. Since only volume preserving variations are 
allowed, 5S has zero integral over S. We assume that there exists a smooth function 
6F/6T, such that for all variations 6T,: 

Dj:F{v,E) ■ 6j: = f ^6J:dS 
is 

where Dy.F is the derivative of F holding v constant 
The Poisson bracket on functions F, G E A is defined by 

where u = cuilv and 

6F _ 6F 
64) 6v J. 

This last term corresponds to the variational derivative of F taken with respect to variations 
of V by potential fiows. 

It has been checked in [17] that the Hamiltonian equation F = {F, H} is equivalent to the 
equations of a liquid drop 

dv ^ , 

— + V^w = -gradp, — = v-n, divt; = 0, p\j: = tk, 
at ot 

where k is the mean curvature of the surface S and r is the surface tension. The Hamiltonian 

is taken to be 

H{v,T) = - [ \\vf dV + T [ dS. 

However this bracket does not define a valid Hamiltonian structure since it is not closed. 
To show that, we will compute the bracket of two specific admissible functionals and show 
that the bracket is not an admissible functional. Let 

F{v,T.) = \l f{\\vf)dV, G{v,T.) = \l g{\\vf)dV, 



where / and g are smooth real functions. Those functionals are admissible and we have 

6F 1 2^ 6G 1 2. 

where Xf (resp. Xg) is the (L^)-orthogonal projection of the vector field (resp. 
^f'dlull^)?; onto the space of divergence free vector fields. 

Proposition 6.1. {F,G} is not an admissible function. 



^^This definition requires us to extend smoothly w in a neighborhood of E. One can check that 5T, is 
independent on the way v is extended and that it is determined up to an additive constant. 
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Proof. We have 

h{v,j:) = {f,g}{v,e) 

= [ cnT\v{XfxX,)dV + l [ {f{\\vf){X,.n)-g{\\vf){Xfn)}dS 

= [ cm\v{XfXX,)dV+ [ {f{\\vf){vXg)-g'{\\vf){vXj)} dV. 

Let's denote the first integral in this expression by Hi and the second one by H2. We have 
1 f r fd^f , \ fd^g 

■ 



+ ■ {D„X, ■ Sv) ^ -£ ■ {D,X, ■ Sv)}dV 



which can be rewritten as 



D„H,-6v = l f I (il.X^■Sv-(^i■Xr]■Sv 



+ {D,Xg ■ Xf) ■ 6v - {D,Xf ■ Xg) ■ Sv^V, 

using the property of symmetry of second Frechet derivative. That is the partial Frechet 
derivative of H2 relative to v admit a gradient. Therefore, this will be the case for H if and 
only if this is true for Hi. We have 

D^Hi ■5v= I cni\{5v) ■ {Xj x Xg) dV 

+ [ curl V ■ {[D^Xf ■ 5v] x Xg) dV + I curl v-{XfX [D^Xg ■ 5v]) dV. 

In this expression, the last two terms are of gradient type because of the symmetry of the 
second Frechet derivative. The first term can be rewritten as 

[ {6v ■ cuT\{Xf X Xg) + diY{6v X {Xf x Xg))} dV, 

which is definitely not of gradient type because of the divergence term. This achieves the 
proof that {F, G } is not an admissible function. □ 

Remark 6.2. The second LMMR bracket presented in Section 5.2 can also be defined for free 
boundary problems with the same difficulties, that is the necessity of a non trivial boundary 
condition in the definition of admissible functionals. Is it possible to define a usable and valid 
bracket for free boundary problems using the method of Soloviev ? 
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